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Abstract. In this paper we consider finite volume hyperbolic manifolds X with non-empty 
totally geodesic boundary. We consider the distribution of the times for the geodesic flow to 
hit the boundary and derive a formula for the moments of the associated random variable in 
terms of the orthospectrum. We show that the the first two moments correspond to two cases 
of known identities for the orthospectrum. We further obtain an explicit formula in terms of 
the trilogarithm functions for the average time for the geodesic flow to hit the boundary in the 
surface case, using the third moment. 



1. Introduction 

Let X be a finite volume hyperbolic manifold with non-empty totally geodesic boundary dX. 
An orthogeodesic for M is a geodesic arc with endpoints perpendicular to dX. These were first 
introduced by Basmajian in pQ in the study of totally geodesic submanifolds. We denote by Ox = 
{cti} the collection of orthogeodesics of X and let k be the length of aj. We note that Ox is 
countable as the elements correspond to a subset of the collection of closed geodesies of the double 
of X along its boundary. We call the set Lx = {h} (with multiplicities) the orthospectrum. 

In PQ, Basmajian derived the following boundary orthospectrum identity; 

(1) Vol(flX) = 2 V n-i ( lo § (coth l - 

where V n (r) is the volume of the ball of radius r in HP. The identity comes from considering the 
universal cover X C H" of X. Then dX is a countable collection of disjoint hyperbolic hyperplanes 
which are the lifts of the boundary components of dX. For each component C of dX, we orthogo- 
nally project each of the other components of dX onto C to obtain a collection of disjoint disks on 
each component C of dX. These disks form an equivariant family of disks that are full measure in 
dX. They descend to a family of disjoint disks in dX of full measure. As each orthogeodesic lifts 
to a perpendicular between two components of dX, each orthogeodesic corresponds to two disks 
(one at each end) in the family of disks in dX and this gives the above identity. 

Using a decomposition of the unit tangent bundle, Bridgeman-Kahn (see [5]) derive the identity; 

(2) Vol(T x (X)) = J2 H n{l) 

ieL x 



) 
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where H n is some smooth function depending only on the dimension n. As Vol(Ti(X)) = Vol(X).V^_i 
where V n —i is the volume of the unit sphere in K", the above identity can also be written as 

Vol(X) = ]T F ™W 

where H n (l) = H n (l)/V n ^. 

In the specific case of dimension two, the function H n is given in terms of the Rogers dilogarithm 
(see [3]). In the papers 0[7] Calegari gives an alternative derivation of the identity in equation [2j 
The motivation for this paper was to connect the above two identities in a natural framework. 
The connection is that they are the first two moments of the Liouville measure. A second motivation 
was to compute the average time it takes to hit the boundary under the geodesic flow. This can 
be put into the same framework and it turns out that consideration of the third moment gives a 
formula for the average time it takes to hit the boundary of M under the geodesic flow in terms 
of the orthospectrum. It is conceivable that higher moments encode other important geometric 
invariants of the manifold. 



2. Moments of Liouville measure 

We let G(HP) be the space of oriented geodesies in EP. By identifying a geodesic with its 
endpoints on the sphere at infinity, the space G(H n ) ~ (S^T 1 x S^T 1 — Diagonal). The Liouville 
measure /U on G(H") is a Mobius invariant measure. In the upper half space model, we identify a 
geodesic with its endpoints (x,y) £ IR™" 1 x Then the Liouville measure /j, has the form 

2dV x dV y 

m-,v) - \ x _ y \2n-2 

where dV x = dx\dx% . . . dx n -i, for x — {x\, X2, ■ ■ ■ , x n -i) € 

If A is a hyperbolic n-manifold with totally geodesic boundary, we identify A, the universal 
cover of A as a subset of W 1 and L C Isom + (H") such that A = A/L. Then G(A) C G(H" ) is the 
set of geodesies intersecting A. We define G(A) = G(A)/L, the space of geodesies in A. Then by 
invariance of Liouville measure \i descends to a measure on G(A) which we also call /i. 

We define the measurable function L : G(A) — > [0, oo] by L(g) := Length(g), where the length is 
taken in A. This is the hitting length function for A. As the limit set Lr has measure zero, almost 
every geodesic hits the boundary of A and therefore for almost every geodesic g, L{g) is finite and 
g is a proper geodesic arc. 

We define the pushforward measure M := L*(/Lt) on the real line. This measure is the distribution 
of lengths of geodesies in A. We define its k-th moment to be 

M k (X) = M(x k ) = [ x k dM=[ L k (g)dfi. 

Jo JG(X) 

In general, the moments of a random variable give a set of measurements that describe distributional 
properties of the random variable such as the average value and variance. Using a decomposition 
of G(A), we show that the moments Mfc(A) have formulae that extend the identities in equations 

fH2i 



The main result of the paper is the following: 
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Theorem 2.1. (Main Theorem) There exists smooth functions F n ^ ■ K + — > R + and constants 
C n > such that if X is a compact hyperbolic n-manifold with totally geodesic boundary dX ^ 0, 
then 

(1) The moment M^(X) satisfies 

M k (X) = J2 F nA l ) 

leL x 

(2) M (X) = C n .Vol(dX) and the identity for M (X) is the identity in equation!!] 

(3) Mi(X) = Vol(Ti(X)) and the identity for Mi(X) is the identity in equation^ 

(4) M 2 (X) = 2Vo\(T 1 (X))A(X) where A{X) is the average time for a vector in T X {X) to hit 
the boundary under geodesic flow. Therefore by the identity for M 2 {X) 

1 



A ^ = 2voimw) E FnAl) = E Gn(l) - 



x 



In the surface case we obtain an explicit formula for the function Gi and hence A(X) in terms of 
poly logarithms. Furthermore, besides compact surfaces obtained as quotients of Fuchsian groups, 
the identity holds more generally for finite area surfaces, which we describe next. 

If S is a finite area surface with totally geodesic boundary dS ^ 0, then the boundary components 
are either closed geodesies or bi-infinite geodesies with cuspidal endpoints. We define a boundary 
cusp of S to be an ideal vertex of dS. We let Cs be the number of boundary cusps of S. Then we 
have the following explicit formula for A(S): 

Theorem 2.2. Let S be a finite area hyperbolic surface with non-empty totally geodesic boundary. 
Then 



47T 2 

F(a) = -12C(3)-— log(l- a ) + 61og 2 (l-a)log(a)~41og(l- a )log 2 (a) 



whe 



-8 log yY^J Li2 ^ + 24L *s(a) + 12l4 3 (l - a), 
Lik(x) is the k th —polylogarithm function, and £ is the Riemann (,— function. 

3. A NATURAL FIBERING 

We have the natural fiber bundle p : T^H") — > G(M n ) such that v is tangent to the oriented 
geodesic p(v). Let fi be the volume measure on Ti(H n ). We parametrize 7i(H n ) as follows: We 
first choose a basepoint b g on each geodesic g (say by taking the point closest to a fixed point p). 
Then a vector v £ Ti(H n ) is given by a triple (x, y, I) G x ]R™ _1 x R where v is tangent to the 

geodesic g with endpoints x, y (from x to y), and I is the signed length along g from the basepoint 
b g . In terms of this parametrization the volume form Q on Ti(H") is 

_ 2dV x dV y dl _ ^ 

where is the Liouville measure (see [ID]). If X is a hyperbolic n-manifold with totally geodesic 
boundary and X the universal cover of X, then the fiber bundle p restricts to Ti(X) to give the 
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equivariant map p : T\{X) — > G(X) which descends to a map p : Ti(X) — > G(X). We have the 
function L : Ti(X) —¥ [0, oo] given by L = L o p. The measure N = L^il was introduced in [5] in 
order to derive the surface case of the identity [2] The measures M, N have a simple relation which 
we now describe. 

For <f) : [0, oo) -jla smooth function with compact support then 

N{4>) = fi(0oL) = / <(>(L(v))dSl v = I I (j)[L(g))dn g dl v 



<f>(L(g)) / dlA df, g = / 4>(L(g))L(g)dfi g = M(x<f>) 

g<£G(X) \Jvep-i{g) ) J geG(X) 

It follows that the measures M, N satisfy dN — xdM. We define the moments of N to be Nk(X) = 
N(x k ). Then 

(3) N k {X) = M k+1 {X). 

Also as iVoPO = Vol(Ti(X)) it follows that Mt(X) = Vo^T^X)). 

Because of the above relation between the measures M, N, the results described in this paper 
can be given in terms of either. For the most part we will give our results in terms of the measure 
M and its moments, but when it is more natural to do so, we will consider the measure N. 

We let A(X) be the average time for a vector in T\(X) to hit the boundary under geodesic flow. 
Then as the time for v and — v sum to L(v), A(X) is half the average of the function L. Then A(X) 
is given by the first moment of measure N and as Ni(X) = M2(X), we have the formula 

1 M 2 (X) 
K ' 2Vol(T 1 (X)) 

proving part 4) of the Main Theorem. 

4. Moments are finite 

Before we derive summation formulae for the moments Mk(X), we need to first show that 
they are finite. The proof that Mq(X), M\(X) are finite will follow from explicit calculation, in 
particular, from the last section Mi(X) = Vol(Ti(X)) which is finite by assumption. By equation 
§ Mfc(X) = N k -i(X), and therefore we need only show that Nk(X) is finite for k > 1. To prove 
this, we show that the measure N on the real line decays exponentially, i.e. there exist positive 
constants a, C such that dN < Ce~ at dt for t large. Then Nk{X) is finite as the measure x k e~ at dt 
is finite. 

We first recall some background on Kleinian groups (see [S] for details). A Kleinian group V is 
a discrete subgroup of the isometries of EP. The limit set Ly = Tx n S^ 1 is the accumulation 
set of an orbit of a point x on the boundary. It is easy to show that Lr is independent of x. The 
convex hull of T, denoted H(T), is the smallest convex set containing all geodesies with endpoints 
in Lr- As H(T) is invariant under T, the convex core is defined to be C(T) = H(T)/T. A Kleinian 
group is convex cocompact if C(T) is compact. Also a group is geometrically finite if N e (C(T)), the 
e -neighborhood of the core, is finite volume. 

Let r be a convex cocompact Kleinian group with N = M n /T and X = H(T)/T its convex core. 
We let <5(r) be the Hausdorff dimension of the limit set Ly- Let gt be the geodesic flow on Ti(iV) 
and define 

B{t) = {v e TipQ I g t (v) g T X (X)} = g t (Ti(X)) n T X (X) 
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The set B(t) is the set of tangent vectors that remain in the convex core under time t flow. We 
now use a standard counting argument on orbits to bound the volume of the set B(t) (see [TO] for 
background) . 

Lemma 4.1. Given T a convex cocompact Kleinian group, then there exists constants A,T such 
that 

Vol(B(t)) < Ae-^- 1 - 5 ^ 1 
for t > T . In particular if Lr ^ SJ^T 1 then Vol(B(t)) is exponentially decaying. 

Proof: Wc take € H(T) and consider its orbits under T. Then we let 

0(r) = {7 e T |rf(0, 7(0)) < r} and N(r) = #0(r). 

By Sullivan (see [E]), there exists constants A, ro such that N(r) < Ae s ^ r for r > ro. We let 
D be the diameter of X. Given a unit tangent vector v, we denote its basepoint by b(v). We let 
B{t) be a lift of B(t) with basepoints within a distance D of 0. For v £ B(t) then by the trangle 
inequality, gt{v) has basepoint b(g t (v)) such that t — D< d(b(g t (v)), 0) <t + D. Therefore b(g t (v)) 
has nearest orbit 7(0) such that b(gt(v)) € B(j(0),D), the ball of radius D around 7(0). Also we 
have that t - 2D < d(-y(0), 0) < t + 2D. Therefore 

U(t)= |J B(j(0),D) 

0(t-2D)"nO(t+2D) 

Let a; be a distance at most D from 0. We want to bound the visual measure of U(t) from x. 
Each B(j(0),D) € U(t) is a distance between t — 3D,t + 3D from x. Let Ti > 3£>, and restrict 
to i > Ti. We radially project each B(j(0),D) e U(t) onto the S(x,t), the sphere of radius £ in 
hyperbolic space. We label the projection P(7(0), D). Then the area of each projection is bounded 
above by a constant C\ > 0. Therefore 

dN(t + 2D) 



Vis x {U{t)) < 



Vol(S(x,t)) 



We have that Vo\(S(x,t)) = S n sinh™ 1 (t) where S n is the volume of the standard Euclidean sphere 
of dimension (n - 1). Thus for t > T u Vo\(S(x,t)) > L n e {n -^ 1 for some L n > 0. 
Therefore for 

for some constant C. 

In order to obtain the bound on Vo\(B(t)) we integrate the visual measure of U(t) over B(t) 
gives 

Vol(B(i)) = Vol(fl(t)) < Vol(X).Ce-« n -^-^ r »* 
for t > T, giving our result for T convex cocompact. □ 

Corollary 4.2. If X is a compact hyperbolic manifold with non-empty totally geodesic boundary 
then the moments Mk{X) are finite for k > 1. 

Proof: Let 

E(t) = {veT 1 (X) I L(v) € [t,t+l)}. 
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Then from the above lemma E(t) C B(t/2) U B(—t/2). Therefore there are constants a, K > such 
that Vol{E{t)) < Ker at for t > 2T a . Therefore 

OO 

N k (X) < Y,{n + l) k .Vol{E{n)) 

n=0 

is finite for k > by comparison with the series ^ n k e~ an . Therefore Mk(X) is finite for k > 1. □ 



5. Decomposition of the space of geodesics 

We let G(H") be the space of oriented geodesics in HI". By identifying a geodesic with its 
endpoints, the Liouville measure \i on G(HP) is given by 

2dV x dV y 



dfi - 



y\ 



2n-2 



If X is a hyperbolic n-manifold with totally geodesic boundary then we let X be the universal 
cover of X as a subset of H" and F C Isom+(H") such that X = X/T. We let G(X) C G(H") be 
the set of geodesics intersecting X. We define G(X) = G(X)/T, the space of geodesics in X. Then 
by invariance of Liouville measure /j, descends to a measure on G{X) which we also call fi. 

The space G(X) has a simple (full measure) decomposition via orthogeodesics. For a an ortho- 
geodesic we define 

F a = {.9 <= G(X) | g is homotopic rel dX to a}. 

The set of orthogeodesics is countable, so we index our orthogeodesics Ox = {&i} and the sets 
Fi = F ai . As the limit set Lp has zero measure, [JFi gives a (full measure) partition of G(X) with 
respect to [i. We note that the sets Fj C SJ^ 1 x S^ 1 are of the form [D\ x D 2 ) U (D 2 x Z?i) where 
Di,D2 are disjoint round disks. Thus G(X) decomposes (up to full measure) into a countable 
collection of elementary pieces indexed by the orthogeodesics. 
Then 



M k (X) = M(x k ) = f x k dM=( L k (g) dv = J2 [ L k (g) 

JO JG(X) \ JF t 



d/j, 



where L : G(X) — > R is the length of a geodesic. 

We now lift Fi to a set F{ C G(H"). In the upper half-space model, we choose two planes P, Q 
with orthogonal distance U = Length(ai). Then Fi lifts to the set of geodesics which intersect 
P and Q. In particular we can take P,Q to have boundary circles centered at with radii l,e li 
respectively. Lifting the function L to Fi we see that it only depends on the endpoints (x, y) and 
the ortholength Zj. Wc denote it by L(x, y, k). 



J Fi J\x\<lJ\v\>e'i \x-y\ 2n 2 



Integrating over the length parameter (in both directions) we obtain 

AL{x,y,li) k dV x dV y 

l\x\<\ JM; 

Therefore we have that 

M k {X) =Y J F n , k {h) 
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and F n< k is given by the integral formula 

4L(x,y,t) k dV x dV y 



F n , k (t) 



\x\<iJ\ v \>e* \x-y\ 2n 2 



This gives the summation formula for the moment 

M k (X) = F n,k( l ) 

and proves part 1) of the Main Theorem. 

By equation^ M t (X) = N {X) = Vol(Ti(A)). Thus letting E 4 C Ti(X) given by E { = p" 1 ^) 
then the identity for M\{X) is given by 

M 1 {X) = Volp^X)) =J2n(Ei). 

i 

In the paper [5] , we show that this is the identity in equation [2j This proves part 3) of the Main 
Theorem. 

6. Zero Moment identity is Basmajian's identity 
As M Q (X) = fi(G(X)) the identity for M (X) is 

i 

We now prove some properties of the Liouville measure \x needed to evaluate both sides of this 
identity and show that it gives Basmajian's identity in equation [T] 

Lemma 6.1. Let /j be the Liouville measure on G(M n ) and P a plane in HP. For U C P, let 
G(U) C G(HP) be the set of geodesies intersecting U transversely. Then the measure /ip on P 
defined by fip(U) — fi(G(U)) is a constant times area measure on P. In particular there is a 
constant K n > depending only on dimension such that /j,p(U) = K n .Vol(U). 

Proof: In the case of n = 2 this is a standard property of fj, (see [5]). In general we see that /zp 
gives a Mobius invariant measure on the hyperbolic plane P and therefore must be a multiple of 
area measure. Thus fJ.p(U) = K n .Vo\(U) where K n only depends on the dimension n. We calculate 
K n in a later section □ 

A Liouville measure preserving map: We let P be a plane in H" and Gp the set of geodesies 
intersecting P transversely. Then Gp = (D- x D + ) U (D + x D ) where D_, D + are disjoint open 
disks in S^* 1 with <9Z?_ = dD + . We let g = (x(g),y(g)) € Gp be the ordered pair of endpoints of g 
and m(g) = gP\P the point of intersection with P. We project the endpoints of g orthogonally onto 
P to obtain the ordered pair of points (p(g), q(g)) G P x P. We note that m(g) is the midpoint of 
the hyperbolic geodesic arc in P joining p(g), q(g)- We now define a map F : Gp —> Gp as follows; 
We first define F : D_ x D + — > Z?_ x D + and extend it to D + x Z?_ by conjugating with the map 
switching endpoints, i.e. if i(x, y) = (y, x) then F(g) = i(F(i(g)) for g e D + x £)_. If g G D- x D+ 
then we define F(g) = h where h is the unique geodesic (in £)_ x D + ) such that m(h) = q{g) and 
q(h) — m(g) (see figure [lj. By construction, F is a involution and if u is an isometry of HP fixing 
P, then F commutes with u. We note that F is not the action of an isometry on the space of 
geodesies. 

Lemma 6.2. The map F : Gp — > Gp preserves Liouville measure. 




Figure 1 . The involution F in the Klein model 
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Proof: We show that F* (p) = fj, by showing the Radon-Nikodym derivative 

Alternately, f(x) is the function such that for any (j) smooth compactly supported function on 
Gp then 

F*{v){<t>) = J iKx)d(F.(ji)) = J <t>{x)y{x)dn{x) 

Therefore 

= = / <l>{F{x))dii{x) = J cf>(x)u(x)d t i(x) 

As F is an involution F(F(x)) = 1, then by the change of variables formula we have 

<j>(x)d[i(x) = [ <j>(F(F(x))d(i(x) = [ cj)(F(x))v(x)d^(x) = [ <j>(x)v(F- 1 (x))v(x)dfi(x) 



Therefore 

v(x)v(F~ x (x)) = 1 or v(x)v{F{x)) = l. 
Similarly we have if u is a hyperbolic isometry then u preserves the Liouville measure and 
= [i. If u also fixes P then f 1 commutes with u and F(u(x)) = u(F(x)). Therefore by the 
change of variables formula again we have 

v({u{x)) = v{x). 

Thus combining the above, if there exists an isometry u fixing P such that u(x) — F(x), then 
v{x) — v(F(x)). But as v(x)v(F(x)) = 1 we then obtain v(x) = 1. 

To find such an isometry, we note that geodesic for x € Gp , and y = F(x) we have the four points 
p(x),m(x) = q(y),q(x) — m(y),p(y) all collinear on P. We choose the plane P* perpendicular to 
P which bisects the hyperbolic interval \p(x),p(y)]. Then refection in P* fixes P and sends x to 
y = F(x). 

Thus v(x) — 1 for all x g Gp and therefore = /x. □ 

Corollary 6.3. //a is an orthogeodesic then 

1(a) 



fi(F a ) = K n .V n -x ( log ( coth 2 

Proof: We consider disjoint planes P, Q with perpendicular distance equal the ortholength 1(a). 
Then /u(-F a ) = fx(S) where 5 C Gp of geodesies which intersect Q. We let Gp = (D_ x D+) U (D+ x 
£>_) where C D + . By the above lemma, fi(S) = /i(F(S)). Let B be the orthogonal projection 
of Q onto P. Then by elementary hyperbolic geometry, B is a ball of radius r = log(coth^i). 
The set F(S) is precisely the set of geodesies in Gp transversely intersecting B which we denote 
by G(B). To see this, we note that if g £ S then F(g) intersects P in B giving F(S) C G(-B). 
Similarly if g G G(B) then F(g) is in 5. Thus, as F is an involution F(S) = G(B). 

Therefore by lemma [O] 

i( a y 



tx(F a ) = fi(S) = KF(S)) = n(G(B)) - F„Vol(B) = i^„K-i I log (^coth g 
□ 

We will now prove the second part of the Main theorem. 
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By the above corollary |6.3| we obtain 



H(G(X)) = V n(F a ) = K n V V n -i ( log ( coth 



1(a) 

/ j f-^y^ a J "li / 4 ' n—l I "-"6 I uuuji 2 
a£Ox a£Ox 

For each boundary component Bi of dX we let B[ in H n be a hyperplane which is a lift of Bi . We 
further take a fundamental domain on B[ for the action of T. We let C be the set of geodesies 
which intersect Di transversely such that the geodesies point into X on B\. Let G = UC then we 
see that G is a lift of G{X) (except for a set of measure zero). To see this note that for almost every 
g G G{X), g is a proper arc from a Bi to a Bj where the orientation of g is pointing into X at Bi 
and out at Bj . Therefore g has lift g' in Cj . Also for i 7^ J if g G Cj (1 Cj then g points inward on 
both Bi and Bj. Therefore C; n Cj = for i 7^ j. Also if <7i, #2 G C, are lifts of the same element 
of G(X) then there is a 7 G T with g^ = 7(31). As Di is a fundamental domain for the action of T 
on B\ then (71 , gi must have endpoints on the boundary of Di which is measure zero. 

Using this we can calculate fi(G(X)). We have 

n(G(X)) = Yu(Ci) 



By the above lemma 6.1 jJb{Ci) — (K n /2)Vol(Bi) where the factor of two comes from Cj containing 
half the geodesies in the set G{Bi) (those pointing into X). Therefore 

MG(A)) = ^i^Vol(^) = ^Vol(dA) 

Thus 

^Vol(dA) = K n Vn-i (log (coth^ 
giving Basmajian's identity. □ 



7. Calculating K n 

To calculate the constant K ni we derive the Lebesgue density of [ip at a point p. Let p be at the 
origin of the Poincare model of H™ and let P be the horizontal plane through p. Let B = _B„_i(p, r) 
be a small (n — 1) -dimensional ball in P about p. Then we have 

2dV X dVy 

where V is volume measure on the unit sphere. If we consider geodesies G(B) then if (x, y) G G(B), 
we have \x — y\ ~ 2. Also if we let 9 X be the angle the ray px from x to p makes with the plane 
P, then the set G(B) X = {y G S™ _1 | (x,?/) G G(B)} is a small ball in the unit sphere. In fact 
the set G{B) is the image of B under stereographic projection from x. Therefore G{B) X is an 
(n — l)-dimensional ellipsoid with the axes of B perpendicular to ray px being approximately 2r 
and other axis approximately 2r sin(9 x ). Therefore 

Vol(G(S x )) ~ 2"- 1 Bin^xJ.Vol^iCr)) 

where B%.(r) is a k-dimensional ball of radius r in the unit sphere. Integrating we get 



dVy \ VoljB^ ir)) 

Ig{b) x \x-y\ 2n - 2 ) x 2- 



KG(B)) = 2.1 , v — dV x ~ V / sSn((?x)dK 
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Thus if A is area measure on P then 

The set St = {x \ 6 X = i] is an (n — 2)-dimensional sphere of radius |cos(t)|. Therefore as 
dV x = d6dV St 

In terms of the Gamma function T we have 

Vol(§* 



2tt^ 

r(*±i) 

Giving 

K - ^ 



2 ™-i r («±i) 
Note: We have K 2 = 1. 

8. Explicit Integral Formulae for F n ,k 

In [3] we derive a formula for L in the surface case. Using this we can write 

k 



F2,k{l) = 



1 



log 



y(y-a)(x-l) 



x(x-a)(y-l) 



-dxdy. 



where a = sech 2 (7/2). 

In [5] we consider the case n > 2 where we derive the explicit formula for L and reduce the 
integral formula to a triple integral via an elementary substitution. Using this we can also reduce 
the integral of F n< k to a triple integral of the form 

= V0l( S ^)VQl( S n-3) f _^^ dr f ^ r ^(mm)\ v 



where 6 = 



1-r 2 ' 

9. Moment Generating function 



The moment generating function of a random variable Y is the function fy(t) = E[e tY \ where 
E is the expected value. We define the moment generating function for measure M 



DO 



(t) = M (e te ) = / e ta; rfAf = / e*' i(9) d/x. 
Jo ./G(X) 

It follows from above that 

cfw= E F »(*>o 

for some function F ra depending only on the dimension n. We similarly can define Cx W = N(e xt ). 
The it follows that the two functions are related by 

C£(0 = | (#(*)) • 
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A simple example, the ideal triangle: We consider the case of X being an ideal triangle. 
In this case, it is more natural to consider measure ./V (in particular M has infinite mass). In ^ 
we show that 

12x 2 
dN = — — g — dx. 
sinh x 

and that 

3(fc + 2)!C(fc + 2) 
N k\X) = 2*=i 

where £ is the Riemann zeta function. In particular the average time to the boundary is 

A(x) Nl{x) 9 cn) 

A(X) 2Vol(T 1 (X)) " 2^ C(3) - 

It follows by integrating that 



cf (*) = y o = 12 (c(2, i - *) + *c(3, i - : 

where £(s, i) is the Hurwitz zeta function 

oo 1 

10. The Surface Case 

For the surface case, the identities in the Main Theorem can be written in terms of polylogarithm 
functions. 

Polylogarithms: The k th polylogarithm function Li^ is defined by the Taylor series 



i ' n K 

i=l 

for \z\ < 1 and by analytic continuation to C. In particular 



Lio(z) = y~, Lii(z) = -log(l - z). 



, Li fc -i(^) . . f* Li fc _ x (t) 

L WJ = givmg Li fc (z) ^ 



Also 



Also the functions Life are related to the Riemann ( function by Life(l) = £(fc). 

Below we describe some properties of the dilogarithm and trilogarithm function. They can all 
be found in the 1991 survey "Structural Properties of Polylogarithms" by L. Lewin (see [5]). 

Dilogarithm: The dilogarithm function Li2(z) is given by 



Lb(*) = -/ l0S(1 " t] dt. 



t 

From the power series representation, it is easy to see that the dilogarithm function satisfies the 
functional equation 

Li 2 (z) +U a (-*) = \^h{z 2 ). 
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Other functional relations of the dilogarithm can be best described by normalizing the dilogarithm 
function. The (extended) Rogers dilogarithm function (see [11) is defined by 

R(x) = Li 2 (x) + i log \x\ log(l -x) x< 1. 

This function arises in calculating hyperbolic volume as the imaginary part of R(z) is the volume 
of the hyperbolic tetrahedron with vertices having cross ratio z. 

Also in terms of the Rogers function, various identities have nice form. Euler's reflection relations 
for the dilogarithm are given by 

R(x) + R(l - x) = R(l) = — 0<x<l 

(4) R(-x) +R(-x- 1 ) = 2R(-1) = -— x>0 

6 

Also Landen's identity is 

(5) R (rzr^) = ~ R ( X ) 0<x<l 
and Abel's functional equation is 

(6) R(x) + R(y) = R(xy) + R + R (v<jLl*l 

\ 1-xy J \ 1 - xy 

In [3] , we showed that the orthospectra of a hyperbolic surface satisfies the following generalized 
orthospectrum identity. 

Theorem 10.1. (Bridgeman, [3]^ Let S be a finite area hyperbolic surface with non-empty totally 
geodesic boundary and Cs boundary cusps. Then 

1 \ 7T 2 (6|x(S)|-C s ) 



ie Ls ^sh 2 ^ 12 



Trilogarithm: 

By definition, the trilogarithm function is given by 



Llz{z) = I ™M dt 



t 

The trilogarithm also satisfies a number of identities. 

(7) Li 3 (z) + Li 3 (-z) = \hi z (z 2 ) 

(8) Li 3 (-z) - Lisi-z- 1 ) = - 1 log 3 (z) - \ log(z) 

6 6 

and 

(9) Li 3 (z) + Li 3 (l - z) - Li 3 (l - z- 1 ) = C(3) + \ log 3 (z) + \ log(z) - \ log 2 (z) log(l - z) 

6 6 2 

If li is an ortholcngth of S, we define 



cosh 2 (£ 
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We will often use the spectrum {ai} instead of {k}. In the paper [3J, we studied the measure N 
and derived the following; 

Theorem 10.2. (Bridgeman, [3]^ There exists a smooth function p : R+ x (0, 1) 

1 AC s x 2 



such that 



Furtht 



where 



dN = ps{x)dx = 



4>(x)p{x, a)dx 



L a {x,y) = ^log 



sinh 2 (x) 



^ p(x, Oi) dx. 



4(/>(L a (x,y))L a (x,y) 
(y-x) 2 

y(y-a){x - 1) 



dxdy 



x(x - a)(y - 1) 

The Moment Identities for Surface: If we apply the above theorem to the function <f>(x) = 1 



To find M k (S), we note that M k (S) = N k -i(S). Therefore by the above, we let <j>(x) 



then we recover the identity in theorem 1 10.1 
get 

M k (S) = N k _t(S) 



„fc-i 



to 



c k 1 p s (x)dx = J2 



We then define 

F k (a) = 
Integrating we have 

Therefore we obtain 



a poo 



'-^fdxdy 

[y - xy 

roa x k+1 dx 



— r^-dxdy 

(y-x) 2 , 

log 



4Cc 



00 x k+1 dx 
sinh 2 (a;) 



a />oo 



y{y-a)(x-l) 
x{x-a)(y-l) 



o sinh (x) 



Ofe-2 / 
L JO 



(fc + l)!C(fc + 1) 
2 k 



(y-x) 2 



-dxdy. 



M k (S) =[J2 Fk 

\l£L s 



cosh 2 (|) / 



(fc + l)!C(fc + l) 

2 k-2 



We note that as £(1) = oo, if S has boundary cusps then Mq(S) is not finite. This corresponds to 
the fact that Mq(S) — Length(dS) which is infinite in the case of boundary cusps. 

Functions Fq , F± are given in terms of simple logarithms and dilogarithms respectively. An 
induction argument shows that F k can be written as a sum of polylogarithm functions of order 
at most k + 1. We will calculate an explicit formula for F<i in terms of trilogarithms in the next 
section. This will give us the formula for the average hitting time A(S) for geodesic flow described 
in theorem 12.21 

11. A Somewhat Brutal Calculation 

We will now obtain an explicit formula for the average hitting time in the surface case in terms 
of sums of polylogarithms evaluated at ortholengths. We let F — F^ given by the above integral 
formula. Then 



F{a) = 



log 



y(y-a)(x-l) 
x{x-a)(y-l) 

(y - x) 2 



-dxdy. 



MOMENTS OF HITTING FUNCTION 



15 



Using Mathematica to calculate the indefinite integral first, gives 7,858 polylogarithm terms which 
then need to be evaluated at the 4 limits to give a final total of approximately 30,000 terms. Also 
the terms must be grouped so that evaluation gives a finite limit. As this seems a daunting task, 
we do the calculation directly using hyperbolic relations to simplify as we go along. The calculation 
is somewhat tedious but the final answer surprisingly short. 

For the reader who would rather skip the long and tedious calculation, evidence for its validity 
is given by figure [2] which is a plot of the difference between the polylogarithm formula for F(a) 
and and its values using numerical integration. As can be seen from the plot, the difference is less 
than 10~ 6 indicating they are the same function. 



6.x LIT 7 



Mi 



- 



/ 




Figure 2. Difference Between Numerical Integration of F and Polylogarithm Formula 



We have that for a £ (0, 1) 



a roc 



F(a) 



Decomposing into cross-ratios, we have 



F(a) 



o Ji 



(log 



y(x-i) 

x(y-l) 



log 



x — a 



(y - xf 



log ( 



x{x-a){y-l) 



{y-xf 



-dydx 



-dydx 



o Ji 



log 



y(x-i) 

x(y-l) 



2-log 



y(x-i) 

x(y-l) 



.log 



y-q 

x — a 



log 



y-q 

x—a 



(y - xf 



-dydx 



Under the mobius transformation m(z) — a/z, we let X = m(x),Y = m(y), then by invariance of 
cross ratios, 

y(x - 1) _ (y - 0)0 - 1) _ (m(y) - m(0)){m(x) - m(l) _ (Y - oo)(X - a) _ X - a 
x(y - 1) ~ (x-Q)(y- 1) ~ (m(x) - m(0)(m(y) - m(l) _ (X - oo)(Y - a) ~ Y -a 
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Thus 



a; 



log 



yQ-i) 

x(y-l) 



(y-x) 2 



-dydx 



-IS 

J oo J a 



1 rO log 



X-a 



Y-a 



v(x-l) 
x(y-l) 



{y -xy 



dYdX 



F(a) = 2. 



a poo 



log 



y-a 

x—a 



a poo 



-dydx + 2. 



log 



log 



x—a 



(y-x) 2 



-dydx 



.log 



y-g 

x—a 



(y - x) 2 



dydx 



(y - x) 2 

We write this as F(a) — 2Ii(a) + 2I 2 (a), where h,I 2 are the above integrals. 

In order to calculate the above integrals we will need the following integral equations. 

(10) f X ^^-dx = log{x)log(l-x/a) + Li 2 (x/a) 

J X (1 

(11) P°g( x ) dx = \ og ( x fi og (i- x / a ) + 2log(x)Li 2 (x/a)-2Li 3 (x/a) 

J X (2 



(12) 



/ 



log(a:) logfr a) ^ = 1 log( ^ )2 lQg(fl) _ log(x)jLz2(a , /a) + Li ^ x / a) 

X 



11.1. Integral 1\. 
Lemma 11.1. 

h(a) = - log(l - a) log 2 (a) + log 3 (l - a) - 4 log 

Proof: We decompose 1\ to obtain 



1 - a 



Li 2 (a) — 6Li 3 



a-1 



h 



log(y — a) 2 — 2 log(y — a) log(a — x) + log(a — x) 2 
iy-x) 2 



dydx = Ji — 2 J 2 + Jz 



Integral J^: 

We have 



Then 



/ 



log(y - a) 2 = log(g/ - of 



y-a 



Also by integral equations above 

/ ^^dy = log (I) lo g (y - a) 2 + 21og(, - a)U 2 (l - |) - 2L, 3 (l - J) 



Thus 



Ji = - 



log(l-a) 3 7T 2 



1 



3 3 ow 3 

Integral J 2 : By parts we have 



+ — log(a)+-log(a) J +log - log(l-a) 2 +21og(l-a)Li 2 1 - -)-2Li 3 1 



log(l - a) log(l - a) log(l - z) 



a — x 



a — x 



As above we have 

r io g (o - x) 



1 — x 



-dx 



■ Iog(a - x) log (L-| j-^l^ 



log(a) log(l-a)+Li 2 
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Also 



J log (« 



log(o -x) 1 2 

dx — — — log(a — x) 

a — x 2 



- x) logfl — x) , 1, , , s , ,, , , , , Sr . 

— - — — -dx = -- log(a - xY logfl - a) + log(a - x)Li 2 

a — x 2 



a- 1 



a - I 



Combining we get 

a log(a-x)(log(l-x)-log(l-a)) 



/ 

Jo 



dx — log(a — x)Li 2 



a- 1 



Thus 



■h = - log(a) log(l - ay - log Li 2 



1 — a I \fl-l 



a- 1 



a- 1 



\og{a)Li 2 ( — - )+Li 3 



Integral J 3 : 

Similarly to J\ we have 

J 3 = [ logfa - xf 
Jo 

Therefore 

J 3 = - log(a - a;) 2 log 



G?y I dx 



1-x 



1 



21og(a — x)Li 2 



a log(a - xf 
1-x 



dx 



a - I 



a - 1 



-2^ 



a — x 



a - I 
a 



a- 1 



1 - a . 

- log(l - a) log(a) 2 + 2 log(a)Li 2 

Thus we have the formula for Ii = J x — 2 J 2 + J 3 giving 
/! = (_l2l(ip)i + i og(a) + 1 i og(a) 3 _ log(a) log(1 _ a) 2 + 2 lQg(1 _ fl)Li2 2iz 3 (1 - A)) 

-2 (- log(a) log(l - a) 2 - log (^) L* 2 (^) + 



log(l - a) log(a) 2 + 2 log(a)Li 2 (^) - 2L* 3 (^) ) 



Using the identities |4j [8j we get 
I 1 (a) = -log(l - a) log 
Simplifying further we also get 



1 -a 



41og Li 2 



1 — a I \ a — 1 



- 6Li 3 



Ii(a) = -log(l-a)log 2 (a)+log 3 (l-a)-41oi ' 



1 - a 



Li 2 (a) - 6Li 3 



a - 1 



a - 1 



□ 
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11.2. Integral la- 
Lemma 11.2. 

2ir 2 1 

J 2 (a) = 2C(3) + — log(l - a) + - log 3 (l - a) - log 2 (l - a) log(a) - log 2 (a) log(l - a) 



-41og(a)Li 2 (a) + 4Li 3 (a) - 2Li 3 
Proof Wo have 



1 -a 



2Li 3 (l - a) 



h = 

/o Ji 

Again we decompose into integrals 



loe 



x(y-l) 



log 

(y-x) 2 



x — a 



a poo 



Jl = 



log (^) log(y - a) 



o Ji 



(y - x) 2 



' <■» iog(i=g)iog(y-o) 

log (^) log(o - x) 
/i (y-a:) 2 
■oo log (2=i) log(a - a;) 



a poo 



Jl 



As before we can integrate to get 

/oo 
log 

and 

J 3 



(y - x ) 2 

log(y - a) 



dydx 

dydx 
dydx 

dydx 

dydx 



1 1 



I 



y-i 

a lo g(i^) log(a-x) 



V - a y 



dy 



1 - x 



dx 



We integrate in J 2 , J 3 the simple factor to get 

\og(y - a) 



Combining J 2 , J 3 we get 

J2 + J3 



(y-x) 



-f 

Jo 



dy I dx = / lo. 



1 — x\ f log(l — x) — log(l — a) i log(l — a) 

1 — x 



a — x 



1-x 



a — x 



dx 



+ 



The Roger's normalized dilogarithm is given by i?(x) = -Li 2 (x) + \ log |x| log(l — x) for x < 1. 
We note that 



R'(x) 



1 / log(l — x) | log |x 

2 I, X" 



+ 



1 - x 
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Thus we let Ri(x) = R(f^) and note that 

1W 2 1 1-x a-x 
Thus 

J 2 + J 3 = 2 j log (^r^) #i(x)dx 
Considering J4 we have similarly 

j4 = j~ log ^ ^ M^f)^ dy = jf" log (i^l) - fl ) - log ^ + log ^ 10 

We then have 

Ji + Ja= / log U - ^ + SV a ; 

We now let #2(2/) = R{a/y), then 



2/ - a 



2 \ y-a 

giving 

Ji + .h = lj™ log tt)^ 

Letting x — a/y we have 



Ji + J 4 = 2 / log 



R'{x)dx = 2log(a)R(a) - 2 [ log(a - x)R , (x)da 

Jo 



a 

Similarly let u = (a — x)/(l — x) then x = (a — u)/(l — u) and 



J 2 + J 3 = 2 / log ( - — - ) R'{u)du = 21og(l - a)i?(a) - 2 / log(o - x)R'{x)dx 
Jo \a — uj J Q 

Giving 

J 2 = 21og(a(l-a))i?(a)-4 / log(a-x)i?'(x)dx = 2 log(a(l-a)i?(a)+2 / log(a-x) ( l0g ^ - ^ + ^ 

Jo Jo V x 1 

We note the formula 
G(x,a) = j ^{a-x)\o g {l-x) dx = 

+ log(a(l - x))L* 2 (l - x) + log (^) U 2 (^) + log (-^) (lA* 

a — x 



Lis 



a(l — x) 
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Thus taking limits we have 

ra log(a - x) log(l - x) 



I 



dx 

o x 



C(3)- — log(a) + -log(a) log 2 (l-a) + log(l-a) \og 2 {a)+log{a{l-a))Li 2 (l-a)+Li 3 (a)-Li 3 {l-a) 
b 2 

From the above 

J 1-x 

- log(x) log lo S [tZ^ ) ~ ^C 1 ~ a ) 1 °g( 1 ~ x ) 1 °g( 2; ) ~ \ ~ a ) lo S 2 (aO 

- log((l - a)a;)Li 2 (a;) - log (j—^J Li ' 2 ( "fTT^" 

, / a — x \ f T ( x — a\ T ( x — a 

- log — — \Li 2 \ - Li 2 



(1 — a)x J \ \ x J \(l — a)x 

T . , , T . ( x — a\ fx — a\ ( x — a 

+ Li 3 {x) + Li 3 + Li 3 [ - - Li 3 [ - r 

\ x J \l — aj \(l — a)x 



Thus taking limits we have 

f a log(a - x) \og{x) 
Jo l ~x 



dx 



ylog(l-a)-ilog 3 (l-a)-log(l-a)log 2 (a)-log((l-a)a)ii 2 (a)+log (rZT^) Li 2 (^A+Li 3 (a)-Li 3 
We now combine to obtain 

h{a) = 2C(3) + ^- log(l log 3 (l - a) - log 2 (l - a) log(a) - log 2 (a) log(l - a) 



-41og(o)Lia(o) + ALi 3 {a) - 2Li 3 ^— — ^ J - 2L» 3 (1 - a) 

□ 

Finally we combine Ii , J 2 to get 

F(a) = 4C(3) + ^ log(l -a) + | log 3 (l - a) - 41og(l - a) log 2 (a) - 2 log 2 (1 - a) log(a) 
-8 log Lz 2 (a) + 8Li 3 (a) - 4L» 3 (1 - a) - 16Li 3 (^) 

Using the identity [9] we get 

47T 2 o n 

F(a) = -12C(3)- — log(l- a ) + 61og 2 (l-a)log(a)-41og(l-a)log 2 ( a ) 

-8 log (t^) Ll2 ^ + 24i *3(a) + 12£i 3 (l - a) 

The function has boundary values F(0) = and F(l) — 12£(3) and is maximized at a = .754493 
with value 17.9804. 

Below is a graph of F. 
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Figure 3. Function F(x) 
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